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Abstract 

The effect of changes in plasma parameters, that are characteristic near or at an L-H transition 
in fusion edge plasmas, on fluctuation correlation lengths are analysed by means of drift- Alfven 
turbulence computations. Scalings by density gradient length, collisionality, plasma beta, and by 
an imposed shear flow are considered. It is found that strongly sheared flows lead to the appearence 
of long-range correlations in electrostatic potential fluctuations parallel and perpendicular to the 
magnetic field. 



This is the preprint version of a manuscript submitted to Plasma Physics and Controlled Fusion. 



1 



I. INTRODUCTION 



The interplay between long-range correlations of turbulent fluctuations, radial electric 
fields and edge bifurcations in 



various experimental studies 



fusion plasmas has received recent interest in the form of 
jl-8]. This interest is motivated by a missing mechanism 
behind the formation of edge transport barriers at the transition from L- to H-mode plasma 
states. The central link between appearence of radial electric field E r and associated sheared 
E r xB flows to suppression of small-scale turbulence, the reduction in turbulent transport, 
and a steepening of the pedestal profile, is generally accepted However, the causal chain 
of mechanisms behind this transport barrier formation is as yet unclear. 

It has been speculated that turbulence generated zonal flows could be able to trigger 
the mean shear flow bifurcation. Long-range correlations in turbulent fluctuations have 
been associated with enhanced zonal flow activity. In L-mode experiments with imposed 
shear flow an increase in correlation length of the fluctuating electrostatic potential |2J and 
density ^| has been found along and across magnetic field lines. 

The influence of single possible players behind the formation of long-range correlations 
can not always easily be determined by experiments, but may be straightforwardly studied 
with numerical simulation. In this work, local drift-Alfven flux-tube turbulence computa- 
tions are applied to analyse correlation statistics for various L-mode parameters in scalings 
that are characteristic for the approach to the H-mode. In particular, scaling effects by 
the background density gradient length, the collisionality, the plasma beta, an imposed 
ExB shear flow and zonal flows on correlation statistics are studied. It is found that only 
strong imposed shear flows are able to generate significant long-range correlations in these 
simulations. 

The work is organized as follows: In Sec. II the numerical model and reference parameters 
are described. In Sec. Ill the evaluation of correlation functions from fluctuating simulated 
quantities is reviewed. In Sec. IV the individual scaling relations are analyzed, followed by 
the conclusions in Sec. V. 
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II. NUMERICAL MODEL: DRIFT-ALFVEN TURBULENCE 



The four-field drift-Alfven fluid model 9] for electromagnetic fusion edge plasma tur- 



bulence is solved numerically using the local flux-tube code ATTEMPT (lOj. The model 
describes the evolution of fluctuations of the electrostatic potential <p, particle density n, 
vector potential A» and parallel ion velocity u: 



dn 1 T 

^ = -Vii J -nV\\u-nlCU) + — K(n) (1) 
at e e 

nrrii dw 

~Bl~dt = V II J + T - /C H ( 2 ) 

T^ + ^^r = — V,| n-0 -77|, J (3) 
at e z n dt en 

du 

nmi ~di = eV H n W 
This set of equations is coupled to the solution of Poisson's and Ampere's equations for 
the vorticity w and the vector potential An: 



W = \7 2 L 4> and fi J = -V 2 ± A\\. (5) 
Operator abbreviations have been introduced as follows: 
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V.ff.V (9) 

A static toroidal equilibrium background magnetic field B is assumed. The model de- 
scribes nonlinear electromagnetic ExB drift motions of electrons of mass m e and ions of 
mass rrii with charges q = ±e. The ion and electron particle densities are equal, n e = = n, 
obeying quasi- neutrality. Ions are cold and electrons have the constant temperature T e , and 
the electron-ion collision frequency is v e . 

A local approximation is applied, where the density gradient is linear and constant in time 
with L~ l = | V In no | with axisymmetric background density n and the density n = n + n 



splitted into a static and a fluctuating part. In the following the tilde on the fluctuating 
density will be avoided for better readability. A partially field-aligned flux-tube coordi- 
nate system {x, y, z) is introduced and the standard drift normalisation is applied, which 



are described in detail in the Appendix of Ref. [10|, where the coordinates (x,y,z) are de- 
noted by (x,rj,a). Parallel derivatives (in z direction) are normalized with respect to the 
parallel connection length L\\ = 2irqRo, perpendicular derivatives with L±, and time scales 
with (c s /L±). The density gradient length then enters via A„ = L±/L n = \d x hxno\- The 
normalized set of equations is: 



OTl 

— = -{(f), n} - X n d y <P - K{<f> -n) + Vj, (J - u) (10) 

^ = -{0,™} + V||J+/C(n) (11) 
8 Aw BT 

+ ^ = -{<M} + V||(n-^)-CJ (12) 

= -{^ u >- v n n ( 13 ) 

In a large aspect ratio circular flux-tube geometry, the Poisson bracket is {f,g} = 
9 x f d y g — d y f d x g, the curvature operator is JC(f) = —ujb [cos(s) d x f + sin(s) d y f], the 
parallel derivative Vy/ = d z f — (3 {A, /}, and the Laplacian becomes = d x + dy. The 
numerical methods using a higher-order Adams-Bashforth / Arakawa scheme are detailed 



in Ref. 



io|. 



Simulations have been performed using reference edge parameters typical of the TEXTOR 
experiment, with major radius Rq = 1.74 m, minor radius a = 0.5 m, electron temperature 
T e = 51.8 eV, magnetic field strength B = 1.0 T, plasma density no = 5.6 • 10 18 m -3 , 
a background density gradient reference scale L± = 3.54 cm, and a parallel connec- 
tion scale L\\ = qR = 465 cm with q = 2.66 and s = (a/q)(dq/dr) = 1. Conver- 



sion to dimensionless model parameters 10| gives a parallel to perpendicular scale ratio 



e = (£||/£±) 2 = 17226, collisionality C = jiL ± /c s v e / 1.96 = 1.0, beta $ = fi n T e /B$ = 1.0, 
mass ratio fx = im e /mi = 4.69, and curvature scale ujb = 2L^/R = 0.046. The numerical 
grid resolution is set to L x x L y x L z = Q4p s x 256p s x 16Ln. 

The magnitude of the dimensionless parameters C, $ and fx in the order of unity is 
typical for many fusion edge plasmas, including larger tokamaks and some stellarator exper- 
iments. The simulations and results are therefore rather generic and not restricted in their 



applicability on a specific tokamak configuration like TEXTOR. 

These nominal values are varied in the simulations to account for changes in pedestal 
parameters related to the approach towards H-mode conditions. The simulations are run 
into a fully developed saturated turbulent state, where time series of density and potential 
fluctuations are recorded at several "probe" position, and are subjected to a correlation 
length analysis. 



III. EVALUATION OF CORRELATION FUNCTIONS 

In this section correlation functions used in the following analysis are reviewed. The 



auto-correlation (AC) function ^auto of a fluctuation signal fit) is defined as [111 ]: 



T 



lauto{r) = 1 t= ° JJ-, • (14) 



E^+ r ) /(*) 

t=0 

T f{tf 

A windowed AC analysis on the computed time series is performed by shifting a slice of the 
data / of size AT by St for every step. Here we use 5t = 0.3 L±/c s and AT = 60 L±/c s . 
The AC function is evaluated in the interval [ti, t\ + AT] with to — 0, t± = 8t, U = % 8t' 

AT+i S T 

E f(t + r)f(t) 

lautoir, U) = (15) 

By evaluation at every time step a time series of the self correlation time TAc{ti) defined by 
Jauto( T Ac(ti),ti) = 0.5 is obtained. The statistical properties of tac(U) are displayed using 
probability density functions (PDF): 

P{tac) = P(h-i < t ac <t b = + dt b ) = — T ^c) (16) 

tb-l<T~AC<~tb 

where N is the length of the fluctuation time series f(t), t b is the position of a bin center, 
with dtb — ( maxr Ac — m i nr Ac)/Nb the width of a bin and N b the number of bins used. P 
gives the probability of finding the auto correlation time r^c i n the fluctuation time series. 

Spatial correlation lengths are analysed by means of the cross-correlation function of two 
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time series f(t) and g(t) are fluctuation time series at two spatial positions: 

^ £(/(* + r) -/) ( S (i) - ») 

73/(T) = »(/) o l9 ) < 17 > 

with 

T T 

f=fi2f and °(f) = T\J2(f-f) 2 - ( 18 ) 

t=0 \ t=0 

To get a measure for the spatial coherence of fluctuations the cross-correlation coefficent 
CC(f,g) = 7 9 /(0) is evaluated as the correlation function 7 s /(t) in the limit r = 0. The 
spatial correlation function Lqc is calculated as the cross-correlation between a fluctuation 
signal / taken at a probe at position lo'f(lo,t) and at a spatially shifted position lj-.f(lj,t), 
with the distance between the probes 51 := \lj — Iq\: 

L cc (5l) = CC(f(l ),f(l + 5l)) (19) 

A statistical description is used, where the data f(t,l) is cut into pieces of length AT, 
as for the auto-correlation time above. The correlation length (Tl9|) is evaluated for a time 
window [iSt ■ AT + i5t]- A time series of half width times Xi(U) results with i = 0,£i = 
S T ,ti = iS T : 

L CC (Sl, U) = CC{f(U, l ), f{U, l + 51)) (20) 

The correlation length A;(tj) is defined as the half width of the correlation function at 
Lcc{^i{ti),ti) = 0.5. Xi(ti) is binned into a histogram and normalised to one, giving a 
probability density function 



P(A,) = P(Z 6 _x < < l b = Z b _! + dl b ) = 1 ^ 5(t-r AC ) (21) 
where Z&_i, are bin centers, cZZf, is the width of the bins. 



N 



IV. SCALINGS OF CORRELATION LENGTHS 



When the L-H transition is approached from an L-mode state, several parameters of the 
edge pedestal are changing that characteristically influence the turbulence and transport. In 
the transition to the H-mode the pedestal density and temperature rise. The collisionality 
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C ~ \ jv e is reduced as the edge temperature grows, the plasma beta (3 = enT/ fioB 2 increases 
and the density gradient length L n becomes smaller. 

Around the L-H transition a mean E x B flow shear layer would develop within the 
edge pedestal region. As turbulence codes to date are unable to self-consistently account 
for realistic H-mode shear flow development, we model this effect by imposing a background 
vorticity on the turbulence. 

The influence of these respective parameter scalings, which model the approach to an H- 
mode state, on fluctuation correlation statistics is analysed in the following. The reference 
"probe" position, at which the time series are recorded, is located in the center of the 
computational domain, corresponding to mid-pedestal radius {L/2) at the torus outboard 
midplane position. Further analyses have been performed for a number of radial "probe" 
position {L/4, 3L/8, 5L/8, 3L/4), which showed very similar results concerning scaling 
relations compared to the radial reference position. Therefore only results for this reference 
position are presented. 

A. Density gradient length scaling 

First, the steepening of the edge density gradient is modelled by varying the density 
gradient length L n while all other parameters remain constant at their nominal L-mode 
levels. A reference simulation is performed with initial gradient length L n0 = 3.54 cm, 
and four simulations with steepened gradient lengths L±/L n = (1.25, 1.5, 1.75, 2.0)Lj_/L n , 
corresponding to physical gradient lengths L n = (2.99,2.49,2.14, 1.87) cm. 

In figure (ED a) global averages (over the whole computational domain except boundary 
dissipation regions) of energetic quantities are shown. The global mean is in addition aver- 
aged over time during the saturated turbulent phase of the simulations, and the temporal 
standard deviations of the fluctuating global quantities are shown as error bars. The values 
are normalised with respect to the reference simulation L n = L n $. 

The zonal flow strength v 2 ZF = ((d x &o) 2 ) x , with $>o(x) = (4>) y , z , is doubled when the 
gradient is steepend corresponding to half the reference gradient length. The zonal flow 
shear = (c^$o) increases by an order of magnitude, similar to the average free energy 
density (\n\). The radial density transport, defined as T x = uvexb,x ~ n(d y (p), increases by 
nearly two orders of magnitude. The enhanced gradient drive is thus found to increase all 
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FIG. 1: Global averages of the fluctuation kinetic energy (Eexb) = ( v | x s) (thin solid 
line), zonal flow energy ((<9 X 0) 2 ) (bold solid line), radial transport (T x ) = (vExB,xn) (dash 
dotted line), density free energy (\n\) (bold dash dotted line), and average zonal flow shear 
(Q) = (dx<p) (bold dashed line). Scalings with varying (a) density gradient length L n , (b) 
collisionality C, (c) normalised plasma beta (3, (d) imposed flow shear Qq — (dvzF/dx)o- 
Temporal standard deviations of the fluctuating global quantities are shown as error bars. 



turbulent activities. 

In figure ([2] a-1) the PDFs of the fluctuation auto-correlation P(tac) are shown for 
electrostatic potential perturbations <fi on the negative axis, and density perturbations n 
on the positive axis. Maxima of the auto-correlation time are found around 4 L^/c s . For 
increasing density drive, secondary smaller peaks emerge around 1 L±/c s . The mean value 
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of density AC is lower by around 0.5 L±/c s compared to potential AC. Mean values of AC 
times are drawn as vertical lines. 

Figure (T3] a-2) shows on the negative and positive axes respectively the PDFs for two 
"typical" nonlinear time scales. 

To estimate time scales for turbulent processes, we first introduce a "convective density 
time scale" r conv which should serve as a measure for the rate of change of the density, given 
by the continuity equation (which in the turbulent state is mainly determined by nonlinear 
convection): d t n ~ v ExB • Vn = njr conv such that T conv = {n/d t n) = {n/\v ExB ■ Vra|). 

As a second time scale, a measure of the local radial ExB velocity is introduced through 
v ExB ~ d y (p = Ps/t ExB , giving r ExB = p s /d y (p in normalised units. 

The PDF with respect to r conv is shown on the right half-space of fig. ([2] a-2): no scaling 
of this time scale with gradient length is observable. The PDF with respect to r ExB is shown 
on the left half-space of fig. ([2ja-2): increasing the gradient drive (corresponding to smaller 
L n ) leads to a maximum of the PDF at smaller t E xb, which therefore can be interpreted as 
enhancing the radial ExB velocity. The correlation time of this maximum corresponds to 
the emerging second maximum in the auto correlation statistics of figure (T5]a-1). 

For the correlation length analysis fluctuations of density n and electrostatic potential <f> 
have been recorded at three further positions, once radially separated by Ax , once poloidally 
(or rather, perpendicular to radius and magnetic field in the drift plane) by Ay, and also 
along the magnetic field line by Az. Correlation length PDFs P(X X ), P{\), P(^z), as 
defined in eq. ( 12T1) . are shown in fig. ([3] a) for increasing density gradient L n . 

The potential correlation lengths, drawn as PDFs on the negative axis of plots ([3] a), 
show a slight rise in correlation lengths radially and poloidally, whereas the correlation 
length along the magnetic field line is slightly reduced when L n is increased. The density 
PDFs (on the positive axis) show an increase in events with short spatial correlations, both 
radially and in parallel, whereas for larger scales the probability is reduced. The density 
PDFs in poloidal direction remain nearly unchanged. 

All spatial correlation length PDF appear to be far from Gaussian, with a steep peak 
at small scales and long tails for larger correlation lengths. The perpendicular scales are 
consistent with dominant turbulent vortex structures of the order of a few p s k±. The net 
effect from an increased density gradient in the constant in general are slightly smaller spatial 
correlation lengths for density fluctuations. 
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FIG. 2: Temporal correlation scalings for (a) L n , (b) C, (c) (3, (d) Qq: (1) auto-correlation 
PDFs for density (positive axis) and potential (negative axis); (2) PDFs of density 
convection time scale r conv (positive axis) and E x B time scale texb (negative axis). 
Mean values of AC times are drawn as vertical lines. 
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TABLE I: Parameters for simulations with increasing collisionality C 
B. Collisionality scaling 

Next, the effect of reduced collisionality on edge turbulence correlations is analyzed. The 
collisionality C oc 1/T e scales inversely with the electron temperature T e . In table ([T]) the 
parameters used for this simulation series are summarised. A collisionality parameter C = 1 
corresponds to a roughly doubled temperature compared to C — 10. 

In fig. (Tj] b) the energetics scaling with collisionality is shown. Except for the zonal flow 
energy (v 2 ZF ) all global energy averages increase linearly with collisionality. For the radial 
density transport T x not only the mean value, but the fluctuation width (drawn as vertical 
deviation bars around the mean) increases. 

In fig. (T5] b) AC times of density and potential fluctuations are shown. Mean values of 
AC times, drawn as vertical lines in fig. ([2]b-l), show no clear scaling with collisionality. 
Only the E x B AC time, on the negative half in fig. (T2]b-2), shows a clear decrease for 
increasing C. 

All time scales are as usual normalized to the drift drift time scale L±/c s . In this series 
of simulations, however, Lj_/c s is varied (according to tabled]) in addition to C. The reverse 
scaling of texb with T e , which is evident in fig. ([2]b-2) in this T e -dependent normalization, 
would also be qualitatively preserved if te x b were plotted in physical units (with the T e 
dependent normalization eliminated). 

Spatial correlation lengths, shown in fig. (I3]b), do not reveal any clear change in the 
correlation of fluctuation signals with collisionality. 
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C. Plasma beta scaling 



In the H-mode the plasma pressure in the pedestal is elevated. For constant magnetic field 
strength, the plasma beta /3 = n 2 T 2 / /Xq/ B 2 rises. This motivates the following simulation 
series, where the magnetic beta is increased, while keeping the collisionality constant. The 
variation of the electron temperature T e , the particle density n, the time and space scales 
for the various simulation runs are listed in table (TIT]) . 

The global energetics are shown in fig. ([T]c). A reduction of zonal flows (drawn as solid 
line) to about a third for the /3 = 10 compared to the reference parameter case — 1), is 
caused by the enhanced Maxwell stress 12]. 

The zonal flow shear, the density free energy, the E x B flow energy as well as the radial 
E x B density transport increase with beta. 

In fig. (0C-1) the AC PDFs P{rAc{ n ))i drawn in the positive half space, show a slight 
decrease of tag = 3 Lj_/c s down to tac = 2 L±/c s for beta rising from 1 to 10. 

Comparison with the E x B and convective timescales on the negative and positive sides 
of fig. ([2] c-2) shows that the lowered AC time for density fluctuations is accompanied by 
an an increasing turbulent E x B drift velocity. The convective time scale shifts to larger 
values for rising magnetic beta. 

The perpendicular spatial correlation length PDFs for the potential, P(X X ) and P{X y ) 
in fig. ([3] c-1,2), show growing tails (negative axis). The perpendicular size of potential 
structures grows with beta. Together with an increasing fluctuation energy this results in 
nearly unchanged mean decorrelation times tac, refered from the mean of P{tac)i drawn 
on the negative half-space of figure (j2]c-l). 

The size of density perturbations increasingly differs from the potential perturbations 
due to the higher non-adiabaticity (via magnetic flutter) of the electrons. The result is an 
average density correlation length of X x ~ X y ~ 4p s for the density, and \ x ~ 7p s and 
X y ~ 6p s for the potential at /3 — 10.0. 

Along magnetic field lines P(X Z ) drawn in fig. ([3] c-3) shows a clear decrease in correlation 
lengths of both density and potential, caused by the enhanced magnetic flutter. 
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TABLE II: Parameters for simulations with increasing magnetic (5 
D. Imposed shear flow scaling 

Further, a simulation series has been performed with the aim of analysing the impact of 
an imposed sheared E x B mean flow on correlations. 

An external electrostatic zonal potential field 00^) is applied through the nonlinear 
advection operators as [(f), f] — > [(f) + 4> , /]: 

o (x) = (l/2)fi (x + L x /2) 2 (22) 
v {x) = n (x + L x /2) (23) 
d x v (x) = n (24) 

The application of 4>q results in a radially increasing E x B mean drift flow with velocity 
v and constant flow shear tt = (0, 250, 500, 1000, 2000). 

In fig. (CO d) the turbulent transport and fluctuation amplitudes show a reduction for all 
levels of an imposed flow shear. The zonal flow amplitude increases strongly for moderate 
Qq = 250 — 500 due to enhanced zonal vorticity coupling Vq ~ (Q)(R e ) of the Reynolds 
stress drive. For larger Q the zonal flows appear strongly reduced, when the Reynolds 
stress R e = d x (f> d v 4> is lowered by the quenched fluctuation amplitudes. 

The fluctuation AC time PDF for the shear flow experiment in fig. (El d-1) shows for 
Qo = 250 — 500 a shift of the maximum to higher AC times (around 3 — 4 L±/c s ) for 
both density and potential fluctuations. For higher values of Qq the maximum of the PDF 
shifts to smaller tac, an d peaks sharply for the density along the positive axis of plot (T5] 
d-1). Vertical lines indicate mean values {^Acju which reflect first the trend to longer self 
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correlation times and then a drop to shorter living structures. 

To properly reflect the lifetime of a turbulent structure, the AC statistics should be 
computed from a fluctuation time series taken in a co-moving frame, or with the mean flow 
velocity subtracted from the field. At a fixed probe position (which has been applied here 
for consistency with experimental measurements) the AC statistics not only maps the eddie 
turnover, but also the background advection of the structure. The drop in the AC PDF in 
fig. ([2] d-2) is therefore partly debted to the faster decay of perturbations, but also to the 
higher convective E x B velocity of convection. 

Characteristic Ex B convection times texb drawn as PDFs in the left half space of fig. ([2] 
d-2) increase for growing Q . This suggests that convection is mainly caused by the fixed 
background E x B mean flow, whereas the convection caused by E x B drift fluctuations 
are suppressed by the mean flow. Along the positive axis the PDFs flatten for increasing 
Q . 

Spatial correlation PDFs are shown in fig. (E]d). For Qq = 250 — 500, the perpendicular 
functions P(X X ), and P(X y ) do not change significantly. Further increasing Qq gives a shift 
of radial correlation lengths to smaller values, for density from around 4 p s to 3 p s . 

The poloidal correlation function P(\ y )(n) on the right half space of fig. ([3] d-2) increases 
to around 3.5 p s for maximal imposed flow strength, and the potential correlation PDF 
along the negative axis show an increase to 6 p s . PDFs in fig. (j3]d-l,2) of perpendicular 
correlation lengths show that the average length for Qq = is larger radially than poloidally 
for the density. For the imposed shear flow of amplitude Qq = 2000 perturbations of density 
as well as of <p are strained out poloidally and radially quenched. 

Fig. ([3] d-3) shows a slight reduction of parallel correlation length PDFs P{\ z ) with 
initially rising flow shear amplitude for both potential and density . For the largest values 
of the imposed flow shear (drawn as black dash dotted line) the parallel correlation length 
of 4> increases from 6 m to about 8 m. The density PDFs P(X z )(n) shows only a damping 
of events with a very low parallel correlation length X z < 2.5 m. 
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V. CONCLUSIONS 



Correlation lengths of density and electrostatic potential fluctuations for conditions rel- 
evant to an L-mode tokamak edge plasma near the L-H transition have been analysed by 
numerical simulation of drift-Alfven turbulence. Five parameter scalings have been inde- 
pendently performed: density gradient length L n , the collisionality C, the plasma beta (3, 
and by imposing a flow shear Q . 

A reduction of L n (corresponding to a pedestal profile steepening) results in slightly 
enhanced perpendicular correlations lengths for 0, and reduced correlation lengths along 
the magnetic field lines. 

Reducing collisionality C (corresponding to a rise in pedestal temperature) did not show 
any clear and significant scaling of all correlation lengths. 

Increasing the plasma beta (3 has different effects on density and potential correlations. 
Perpendicular correlation lengths for <fi increase, whereas the density correlation length PDF 
is shifted to smaller spatial scales. Parallel correlation lengths are reduced for both n and 
0. 

Externally imposing a flow shear Q was found to significantly enhance poloidal and 
parallel correlations lengths of only for very strong shearing rates. Density correlation 
lengths are increased poloidally but are reduced along magnetic field lines. Radially a 
reduction of correlation amplitudes of n and has been found. 

It can be concluded from these parameter scaling simulations, that experimentally ob- 
served long-range correlations near or at the transition to H-mode states are likely caused 
by the straining effect of a strongly sheared flow on the turbulence. The zonal flows are 
amplified strongly only for moderate imposed mean flow shear, whereas long-range correla- 
tions appear only for strong external shearing. All other plasma parameters scalings, which 
appear towards a transition, have either a weak or reducing influence on correlation lengths. 
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